Localization and Spreading of Diseases in Complex Networks 
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Using the SIS model on unweighted and weighted networks, we consider the disease localization 
phenomenon. In contrast to the well-recognized point of view that diseases infect a finite fraction 
of vertices right above the epidemic threshold, we show that diseases can be localized on a finite 
number of vertices, where hubs and edges with large weights are centers of localization. Our results 
follow from the analysis of standard models of networks and empirical data for real-world networks. 
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Survey of infectious diseases reveals that before an out- 
break, often, if not typically, a disease is localized within 
a small group of individuals. Changes in environmental 
conditions or increase in the frequency of external con- 
tacts result in an epidemic outbreak. In the present pa- 
per we propose an approach that enables us to describe 
quantitatively this important localization-delocalization 
phenomenon. Our approach is based on the SIS model 
[H, Q of spreading of diseases in weighted and unweighted 
networks, where the weights of edges encode frequency of 
contacts between vertices. It is widely accepted that in 
uncorrelated networks the epidemic threshold A c of the 
infection rate A is Xmf = (q) /(q 2 ), where (q) and (q 2 ) are 
the first and second moments of the degree distribution 
[2| . So in networks with a finite (q 2 ) the threshold should 
be non-zero, while it is zero if (q 2 ) diverges. One should 
stress however that all these well-known results were ob- 
tained only within a mean-field theory, actually within 
an annealed network approximation in which a random 
network is substituted for its fully connected weighted 
counterpart @]. Contrastingly, one can show exactly for 
an arbitrary graph that A c is actually determined by the 
largest eigenvalue Ai of the adjacency matrix Aij of the 
graph, and A c = 1/Ai |3-|§|. For uncorrelated networks, 
it was found that Ai is determined by the maximum de- 
gree q ma x, Ai oc y/q m ax 131 ■ Then, if in the infinite size 
limit, q max tends to infinity, as, e.g., in the Erdos-Renyi 
graphs, this leads to an amazing conclusion that the epi- 
demic threshold is absent even in (infinite) networks with 
a finite (q 2 ) in contrast to the mean- field result. 

In the present paper we develop a spectral approach to 
the SIS model on complex networks. We show that the 
contradiction between the mean-field approximation and 
the exact result can be resolved if we take into account 
localization of diseases. It turns out that, in contrast to 
the mean field theory, in which a finite fraction of vertices 
are infected at A > A c , there are actually two scenarios 
of the spreading of diseases. If Ai corresponds to a lo- 
calized eigenstate, then, at A right above A c = 1/Ai, 
disease is mainly localized on a finite number of vertices, 
i.e., the fraction of infected vertices is negligibly small 



in large networks. With further increase of A, disease 
gradually infects more and more vertices until it will in- 
fect a finite fraction of vertices. In the second scenario, 
Ai corresponds to a delocalized state. Then already at 
AAi — 1-Cl, disease infects a finite fraction of vertices. 
Analysing network models and real-world networks, we 
show that hubs, edges with large weights, and other dense 
subgraphs can be centers of localization. 

We consider the standard SIS model of disease spread- 
ing in a complex network of size N having the adjacency 
matrix with arbitrary entries Ay > 0. Infected vertices 
become susceptible with unit rate, and each susceptible 
vertex becomes infected by its infective neighbor with the 
infection rate A. The probability Pi(t) that vertex i is in- 
fected at time t is described by the evolution equation 

^ = - Pl (t) + A[l - Pl (t)} J2 A ijPj (t). (1) 

3=1 

In the steady state, at t — > oo, the probability Pi = Pi(oo) 
is determined by a non-linear equation, 



1 + 



(2) 



which has a non-zero solution pi > if A is larger than a 
so-called epidemic threshold A c . In this case, the preva- 
lence p = 52j=i Pi/N is also non-zero. This transition 
contrasts to a "thermodynamic phase transition" , since 
the critical feature at A c is present even for finite nets. 

Spectral approach. — To solve the SIS model, we use 
the spectral properties of the adjacency matrix A. The 
eigenvalues A and the corresponding eigenvectors / with 
components fa are solutions of the equation A/ = Af. 
Since the matrix A is real and symmetric, its N eigenvec- 
tors /(A) (A max = Ai > A 2 > . . . Ajv) form a complete 
orthonormal basis. According to the Perron-Frobenius 
theorem, the largest eigenvalue Ai and the corresponding 
principal eigenvector /(Ai) of a real nonnegative sym- 
metric matrix are nonnegative Q. The probabilities Pi 
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can be written as a linear superposition, 



= 5>(A)/i(A). 



(3) 



The coefficients c(A) are projections of the vector p on 
/(A). Substituting Eq. © into Eq. fl2J, we obtain 



c(A) 
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/i(A)/i(A') 
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(4) 



In order to find the epidemic threshold A c and p(X) near 
A c , it is enough to take into account only the princi- 
pal eigenvector /(Ai) in Eqs. ((3]) and (j4}, i.e., pi ~ 
c(Ai)/j(Ai). Solving Eq. (gj with respect to c(Ai) in 
the leading order in r = AAi - 1 < 1, wc find A c = 1/Ai 
and p w air, where the coefficient ai is 



A r 



N 



«i=E^( a o/[^E/'( a i)] 



(5) 



This expression is asymptotically exact if there is a gap 
between Ai and A 2 (see also Ref. [loj]). Thus, at r <C 1, p 
is determined by the principal eigenvector. Contributions 
of other eigenvectors are of the order of t 2 . Considering 
two largest eigenvalues in Eq. Ai and A 2 , and their 
eigenvectors, we obtain p(X) ~ a%T + o^r 2 and so on. 

The usual point of view is that a\ is of the order of 
0(1), and so a finite fraction of vertices is infected right 
above A c . To learn if another behavior is possible, we 
study whether Ai corresponds to a localized or delocal- 
ized state. We introduce the inverse participation ratio 



N 



/Pi2(A) = £#(A). 



(6) 



If, in the limit N -> oo, IPR(A) is of the order of 0(1), 
then the eigenvector /(A) is localized. If IPR(A) — > 
then this state is delocalized. For a localized /(A) the 
components /i(A) are of the order of O(l) only at few 
vertices. For a delocalized /(A) we usually have fi(A) ~ 
0(1/ VW) < 1. From Eq. © it follows that if the prin- 
cipal eigenvector f(Ai) is localized, then a% ~ 0(1/N) 
and so p ps ct\T ~ 0(l/N). In this case, above A c disease 
is localized on a finite number Np of vertices. If /(Ai) 
is delocalized, then p is of the order of O(l), and disease 
infects a finite fraction of vertices right above A c . These 
two contrasting scenarios are shown in Fig. [T]for the SIS 
model on the karate-club network ll| and the weighted 
collaboration networks of scientists posting preprints on 
the astrophysics archive at arXiv.org, 1995-1999, and the 
condensed matter archive at January 1, 1995 - March 31, 
2005 The astro- ph and karate-club nets have delo- 
calized principal eigenstates while the cond-mat-2005 net 
has a localized principal eigenstate. Numerical solution 
of Eq. @ gives a± = 1.8 xlO -3 for the astro-ph net and 
ai = 1.5xl0~ 4 for the cond-mat-2005 net. 




FIG. 1. Prevalence p versus the infection rate A in real net- 
works, (a) astro-phys (upper line) and cond-mat-2005 (lower 
line) weighted networks [from Eq. (b) Karate-club net- 

work. The lower curve accounts only the eigenstate Ai in 
Eq. Q. Adding the eigenstate As, we obtain the higher curve 
and so on. The most upper curve is the exact prevalence. 



There is an iterative method to find Ai and IPR(Ai), 
Ai = lim Ai(n) = lim (g (n) !g (n) )/|g (n) | 2 , (7) 
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IPR(A 1 )= lim ^(.g| n) ) 4 /l9 ( " )|4 



(8) 



where = Ag^ n ' and is a positive vector. Ai(n) 

is the lower bound of Ai. For unweighted networks, i.e., 
Ay = 0, 1, choosing g(°'=l, we obtain that the first iter- 



ation n — 1 gives gf^ = qi and 

. i v-^ a a (q)®' 2 "' 

Al(1) = JTfl^N 1^ qiAi i q i = Amf ' 
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2N /(g) and r is the Pearson coefficient 13 1 
14| ) . Eq. shows the influence of degree- 



where Amf= 
(see also Ref. 

degree correlations on Ai. The first iteration also gives 
the mean-field result I PR = (q*) / [N (q 2 ) 2 } - 0(l/N). 
Already a few iterations give good approximations for 
Ai and IPR if the principal eigenstate is delocalized but 
more iterations are needed if this eigenstate is localized. 

Bethe lattice. — What structural elements of a network 
are responsible for the localization of the principal eigen- 
state? Let us consider a Bethe lattice (see Fig. [2]) . The 
adjacency matrix of an unweighted regular Bethe lattice 
with vertices of degree k has the largest eigenvalue Ai=fc. 
The principal eigenvector /i(Ai)=iV _1 / 2 is delocalized. 
Other eigenstates are not important for us. Let us intro- 
duce a hub of degree q>k connected to the neighbors by 
edges with a weight w > 1 [see Fig. EJb)]. Other edges 
have weight 1. Solving the equation Af = Af , we find 



Ai = qw 2 ! / ' V 1 qw 2 — B, 

IPR = f*(A 1 )[l + qw 4 /(a' l -B)}, 

f (A 1 ) = [(qw 2 /2-B)/(qw 2 -B)}^. 



(10) 

(11) 
(12) 
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FIG. 2. (a) Regular Bethe lattice with degree k = 3. (b) 
Bethe lattice with one hub of degree q > k. This hub is 
connected to neighbors by edges having the same weight w > 
1 (red lines), (c) Bethe lattice with two vertices of degrees q\ 
and q2 connected by an edge with a weight w > 1 (red line). 
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Here B = k — 1 is the branching coefficient of the graph, 
a = [qw 2 — B] 1 / 2 , and /o(Ai) is the component of the 
principal eigenvector f(Ai) at the hub. /i(Ai) depends 
only on distance n from the hub to vertex i, i.e., /i(Ai) = 
/n(Ai), and exponentially decays with n > 1, 



/„(Ai) = Wo(Ai)/V- 



(13) 



This exponential decay leads to a finite IPR, so this 
eigenstate is localized. If qw 2 ^> B, then IPR — > (1 + 
1/q) /4. The second eigenstate with A2 = k is delocalized. 
The criterion Ai > k, leads to the condition that q must 
be larger than q ioc = (B 2 + B)/w 2 . For the SIS model 
on this Bethe lattice, Eqs. J3]) and (TT31 give oc 1/N. 

Now we consider the Bethe lattice with two hubs of 
degrees qi and q 2 connected by an edge with a weight 
w > 1 [see Fig. [UJc)]. Other edges have weight 1. We 
find that two eigenstates with eigenvalues Ai and A2 can 
appear above the delocalized eigenstate with A3 = fc, 



Ai( 2 ) = a± + B/a±, 



4 = 2«i- 



(Qi+Q2+vr)-4QiQ 2 



1/2 



*i(4+Qi) + ^(4+Q 2 ) = 4 - B, 
IPR(A m ) = ^i(a 4 ± +Q 1 )+^ 2 (a 4 ± +Q 2 )]/(ai-B). (14) 

The signs ± correspond to Ai and A 2 , respectively, and 
Qi(2) = qi(2) — B— 1. The components fi(Ai( 2 )) decrease 
with increasing distance n from vertex i to the hubs as 
^1(2) / a ± if vertex i belongs to branches connected to 
hub 1 or hub 2, respectively. ^1(2) are the components 
of /(Ai) at hubs 1 and 2, respectively Their ratio is 
^2/^1 = (4 — Qi)/ ( wa ±)- The criterion for localization 
is given by the condition Ai,A 2 > k. If q% = qi and 
w ^> 1, IPR(Ai) reaches the maximum value 0.5 that 
means localization on two hubs, ^1 = ^2 — > I/a/2- 

Scale-free networks. — To study appearance and prop- 
erties of localized eigenstates in uncorrelated complex 
networks, we use the static model [l5| that generates 
unweighted scale-free networks with degree distribution 
P(q) cx Cq~~t at q > 1. Using software OCTAVE, for 
each realization of a random network of size with the 
mean degree (q) and 7 = 4, we calculated eigenvalues, 



FIG. 3. (a) Prevalence p versus A in a scale- free network of 10 5 
vertices generated by the static model with 7 = 4, (q) — 10. 
The lowest curve accounts for only the principal eigenstate 
in Eq. Q, the next one accounts Ai and A2, and so on. (b) 
Zoom of the prevalence at A near A c = 1/Ai. 



eigenvectors, and IPR(A) of the adjacency matrix. In 
networks of size N = 10 5 , we found that several (typi- 
cally, from one to three for different realizations) eigen- 
states appear above the upper delocalized eigenstate. 
These states are localized at hubs and their properties 
are described well by Eqs. pU)) -([T3 ^) with w = 1 if the 
branching coefficient B in these equations is replaced by 
the averaged branching coefficient B — (q 2 ) / (q) — l. In 
these scale-free graphs, the upper delocalized eigenstate 
Ad is slightly above the mean- field value Amf = (q 2 ) / (q)- 
The maximum degree q m ax fluctuates from realization to 
realization. In realizations with a small q m ax, the prin- 
cipal eigenvector is delocalized and Ai = A^. Therefore, 
the criterion for localization of the principal eigenstate 
at a vertex with degree q ma x can be written as 



Ai 



B > Aw. 



(15) 



The equality here gives the threshold degree qi oc . For 
N = 10 5 , (q) = 10, and 7 = 4, our numerical calcula- 
tions give {q 2 }/(q) ~ 14.1 and A d « 15.1. According to 
Eq. (fT5]) . a localized state appears above A^ if q max is 
larger than qi oc w 214. Since the averaged value of q max 
depends on N, at small A^ the probability to generate a 
graph with q max > qi oc is small JTfjj]. Only sufficiently 
large graphs can have a localized principal eigenstate. 

Fig. [3] represents results of our numerical solution of 
Eq. (|4|) for the SIS model on a scale-free network of size 
N = 10 5 with the principal eigenvector localized at hub 
with q = 323. Equations ([TP ]) -([T3" 1) and © give Ai = 
18.35, IPR = 0.23, and a x ~ 1.4xl0~ 3 . These values 
agree well with the measured values Ai = 18.47, IPR — 
0.21, and a± ~ 1.7xl0 -3 . The eigenvector with A2 is 
also localized at another hub with q = 254. The third 
eigenvalue A3 corresponds to a delocalized eigenvector. 
The first two eigenstates allow to describe p(X) close to 
A c =l/Ai. Accounting for the delocalized eigenstate A3 
gives better results in a broader range of A (see Fig. 
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TABLE I. Characteristics of real-world networks. is size, 7 is the degree distribution exponent, q ma x is the maximum 
degree, qi oc is the localization threshold found from Eq. (|15p . Ai is the largest eigenvalue and Ai(l) is its lower bound, Eq. @, 
respectively. D and A stand for assortative and disassortative mixing. Two last columns represent weighted networks. 



Network 


N 


7 


Qma x 


Hoc 


ii 2 )Ki) 


mixing 


Ai 


Ai(l) 


JPii(Ai) 


Ai 
weif 


IPR{k{) 
dited 


cond-mat 2005 [12] 


40421 


3.0 


278 


2604 


27.35 


A 


51.29 


35.205 


0.0081 


47.63 


0.3415 


hep-th [12] 


8361 




50 


521 


8.687 


A 


23 


10.632 


0.0417 


40.52 


0.3531 


astro-ph [12] 


16706 




360 


5415 


44.92 


A 


73.89 


56.287 


0.005 


33.7575 


0.0525 


power grid [17] 


4941 


exponential 


19 


53 


3.87 




7.483 


3.9 


0.041 






fp5 [18] 


27985 


2.2 


2942 


38610 


211.0 




197.03 


176.3 


0.0035 






CAIDA (router-internet) [19] 192244 


2.7 


1071 


11947 


37.89 




109.5 


42.9 


0.010 






karate club fill 


34 




17 


37 


7.77 


D 


6.72 


6.01 


0.073 







Real networks. — The largest eigenvalue Aj., IPR(Ai), 
and other parameters of a few weighted and unweighted 
real- world networks are given in Table [IJ Note first that in 
all these unweighted real networks the inverse participa- 
tion ratio IPR(Ai) is small that evidences a delocalized 
Aj. Second, in unweighted networks, Ai differs strongly 
from the mean- field value A MF = (q 2 )/(q). Ax < A MF 
in networks with disassortative mixing and Ai > Amf in 
assortative networks. Thus, degree-degree correlations 
strongly influence Ai. A similar observation has been 
made in Q. Qualitatively, this result agrees with the 
lower bound of Ai determined by Eq. ©. So we suggest 
that assortative mixing decreases the epidemic thresh- 
old A c = 1/Ai compared to the mean-field value 1/Aa/f 
while disassortative mixing increases A c . 

Table Q] shows that in contrast to the unweighted hep- 
th and cond-mat-2005 networks, their weighted versions 
have a localized principal eigenvector with a large IPR. 
Localization occurs at vertices linked by edges with large 
weights. In the cond-mat-2005 network, localization oc- 
curs at vertices of degrees 37 and 28 connected by an 
edge with weight 34.3 that is much larger than the av- 
eraged weight w — 0.51. Components of the principal 
eigenvectors decrease exponentially from the center of 
localization in agreement with the solution Eq. (|14[) . In 
the hep-th network, this strong edge has weight 34 larger 
than w — 0.97 and connects vertices of degrees 34 and 
33. In unweighted networks localization does not occur 
because the localization threshold qi oc obtained from the 
criterion Eq. (fT5|) exceeds q ma x- The prevalence p(A) in 
two weighted networks is represented in Fig.rTfa). 

In conclusion, based on a spectral approach to the SIS 
model, we showed that if the principal eigenvector of the 
adjacency matrix of a network is localized, then at the 
infection rate A right above the threshold 1/Ai, disease 
is mainly localized on a finite number of vertices. Impor- 
tantly, a strict epidemic threshold in this case is actually 
absent, and a real epidemic affecting a finite fraction of 
vertices occurs after a smooth crossover, at higher val- 
ues of A. On the other hand, if the principal eigenvector 



is delocalized, the epidemic occurs in the whole region 
above A c = 1/ Ai . We suggest that further investigations 
of real-world networks will give many new examples of 
the disease localization-delocalization phenomena. 
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